Spin Liquid States at the vicinity of metal-insulator transition 
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We study in this letter quantum spin liquid states (QSLs) at the vicinity of metal-insulator 
transition. Assuming that the low energy excitations in the QSLs are labeled by "spinon" occupation 
numbers with the same Fermi surface structure as in the corresponding metal (Fermi-liquid) side, 
we propose a phenomenological Landau-like low energy theory for the QSLs and show that the 
usual ?7(1) QSLs is a representative member of this class of spin liquids. Based on our effective 
low energy theory, an alternative picture to the Brinkman-Rice picture of Mott metal-insulator 
transition is proposed. The charge, spin and thermal responses of QSLs are discussed under such a 
phenomenology. 
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Quantum spin liquid states (QSLs) in dimensions d > 
1 has been a long sought dream in condensed matter 
physics which has not been confirmed in realistic mate- 
rials until rather recently. [l| These states are electronic 
Mott insulators that are not magnetically ordered down 
to the lowest temperature due to strong quantum me- 
chanical fluctuations of spins and/or frustrated interac- 
tion. Various exotic properties have been proposed to 
exist in QSLs. For instance, charge neutral and spin- 
i mobile objects, spinous, were proposed to emerge in 
such electronic states at low temperature accompanied 
by different kinds of (emerging) gauge fields. The spinous 
may be gapped or gapless and may obey either boson or 
fermion statistics. These new particles and gauge fields 
which characterize low energy behaviors of the system do 
not appear in the parent Hamiltonian and "emerge" as a 
result of strong correlation. 

In the past few years, several experimental candi- 
dates for QSLs have been discovered that support the 
existence of gapless fermionic spinon excitations. The 
best studied example is a family of organic compounds 
K-(ET)2Cu2(CN)3 (ET)[3] and Pd(dmit)2(EtMe3Sb) 
(dmit salts) [3|. Both materials are Mott insulators in 
proximity to the metal-insulator transition because they 
become superconductor (ET) or metal (dmit) under 
modest pressure. Despite the large magnetic exchange 
J « 250 K observed in these systems, there is no experi- 
mental indication of long range magnetic ordering down 
to temperature ~ 30 mK. Linear temperature depen- 
dence of the specific heat and Pauli-like spin suscepti- 
bility were found in both materials at low temperature 
suggesting that the low energy excitations are spin-1/2 
fermions with a Fermi surface. 0, Q This Fermi liquid- 
like behavior is further supported by their Wilson ra- 
tios which are close to one. The thermal conductiv- 
ity experiments on the ET salts found a large contri- 
bution to K beside phonons with k/T much reduced be- 
low 0.3 K,@ while k/T approaches to a constant down 



to the lowest temperature in dmit salts. [3| All these ex- 
perimental observations point to the scenario that the 
low lying excitations in these Mott insulators are mo- 
bile fermionic particles (spinous) that form a Fermi sur- 
face like their parent electrons. Besides ET and dmit- 
salt, the Kagome compound ZnCu3(OH)6Cl2, the three 
dimensional hyper-Kagome material Na4lr308 and the 
newly discovered triangular compound Ba3CuSb209 are 
also considered to be candidates for QSLs with gapless 
excitations, [sl-fiol 

Several experiments were proposed to probe mobile 
spinous. For example, giant magnetoresistance like 
experiment was designed to measure mobile spinous 
through oscillatory coupling between two ferromagnets 
via a quantum spin liquid spacer. [11] The thermal Hall 
effect in insulating quantum magnets was proposed as an 
example of thermal transport of spinons, where different 
responses were used to distinguish between magnon- and 
spinon- transports. [l^ The spinon life time and mass can 
be measured through sound attenuation experiment. [l^ 
Despite all these proposals, a generic method to compare 
theoretical prediction of QSLs to experimental data is 
still missing at the phenomenological level. 

The purpose of this letter is to build a general phe- 
nomenological theory for spin liquids with (large) Fermi 
surfaces. We observe that like Fermi liquid states, the low 
energy excitations in these QSLs are assumed to be la- 
beled by the same occupation numbers as free fermions, 
the difference between Fermi liquids and QSLs is that 
the excitations in Fermi liquid are quasi-particles that 
carry both charges and spins, whereas the excitations in 
QSLs carry only spins. In particular, DC charge trans- 
port exists in Fermi liquid states, but is absent in QSLs 
(insulators) . 

The fact that the low energy excitations in these 
QSLs are labeled by the same occupation numbers as 
free fermions suggests that the excitation energy AE = 
E — Eq for these states are also given by a Landau-type 
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expression 
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where = — ^ is the (single) spinon energy mea- 
sured from the chemical potential ^, m* is the spinon 
effective mass and a and a' are spin indices. Sn pa — 
n-pa — Wpg. measures the departure of the spinon distri- 
bution function from the ground state distribution Up = 

^(~^p)- /pp' is the interaction energy between excited 
spinous. A spherical, rotational invariant Fermi surface 
is assumed here for simplicity. In this case we may write 
/p^/ in terms of spin symmetric and spin antisymmet- 
ric components f^^', = /pp/^o-o-' + /pp'O-a'. For isotropic 
systems, /pp?^ depends only on the angle 9 between p 



and p' and we can expand /pp?'' = Yll^a fi^'^'Pi{cos9) at 



s(a) - 



3D and f^^' = Y.Zq fi cos{W) at 2D, where Pi's are 
Legendre polynomials. The Landau parameters, defined 

by 

provide a dimensionless measures of the strengths of the 
interactions between spinous on the Fermi surface, where 
A^(0) is the Fermi surface density of states. The low tem- 
perature properties of the QSLs are completely deter- 
mined by the spinon mass m* and the interaction /pp, 

(or F^''-'^^) as in Fermi liquid theory. 

Notice that the energy functional AE for our QSLs is 
so far identical to that for Fermi liquids. We propose 
and shall demonstrate in the following that the QSLs 
distinguish themselves from Fermi liquids by having a 
strong constraint on the Landau parameters F^. 

We start with the observation that the charge current 
J carried by quasi-particles in Fermi liquid theory (and 
in QSLs) is given by 



J^4(l + f)j(°), 

m* a 



(2a) 



where J^^^ is the charge current carried by the corre- 
sponding non-interacting fermions and d is the dimen- 
sion. For translational invariant systems, the charge cur- 
rent carried by quasi-particles is not renormalized and 
^ = 1 -|- -^[l^. However, this is in general not valid for 
electrons in crystals where Galilean invariance is lost. In 
this case ^ 1 + ^^'^ the charge current carries by 
quasi-particles is renormalized by quasi-particle interac- 
tion. On the other hand, the thermal current Jg is only 
renormalized by the effective mass in Fermi liquid theory. 



T _ ™ t(0) 

m* ^ 



(2b) 



where Jg'' is the corresponding thermal current carried 
by non-interacting electrons. (See supplementary mate- 
rials for details.) Thus, in the special case 1 -I- /d — J> 



while ^ remaining finite, J and Jg suggesting 
that the electronic system is in a special state where spin- 
1/2 quasi-particles do not carry charge due to interaction 
but they still carry entropy (i.e. electric insulating but 
thermal conducting). This is exactly what we expect for 
spinous in QSLs. We note that it is crucial that F^ is 
independent of ^ for this mechanism to work. 

To put our argument in a more quantitative frame- 
work we study the electromagnetic responses of a Fermi 
liquid system with l + F^/d^ 0. The charge and (trans- 
verse) current response functions are given by the stan- 
dard Fermi liquid forms [l3. [lil [l7[ 



x<i(q,t^) 
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and 
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where xodil,^) and xot{ci,^) are the density-density 
and (transverse) current-current response functions for 
a Fermi gas with effective mass m* but without Lan- 
dau interactions, respectively. The longitudinal current- 
current response function xi is related to Xd through 
X<i(q,'^) = (<7V'^^)X;(q7w) and the ac conductivity cr;(t) 
is related to xi{t) by cr;(t)(q,a;) = e'^xi{t){(l,^)/i^, where 
q = \q\. In the singular limit l + Ff /d 0, it is clear that 
higher order q, w-dependent terms should be included in 
the Landau interaction to obtain finite results. Expand- 
ing at small q and lu, we obtain 



l + Ffiq,u)/d 
N{0) 



(5) 



where qt ^ Vx and ^ V are associated the transverse 
(curl) and longitudinal (gradient) parts of the small q 
expansion, a = in the QSLs. Putting this into the 
charge response function Xd, we find that to ensure that 
the system is in an incompressible (insulator) state, we 
must have 7; = 0. The other possibility _Fq 00 implies 
complete vanishing of charge responses in the insulating 
state. 

With this parametrization we obtain for the ac con- 
ductivity at small uj, 
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where (7q{uj) = e2xot(0, i:j)/(ii:j) = e^xo^Oi '^)/(*'^)- The 
last equality is valid as long as F^ is finite. Eq. ([6]) 
was first obtained in the U(\) gauge theory approach 
to spin liquid in a slightly different form [3] and predicts 



power law conductivity Re[(T(a;)] 



3.33 



(w2) (at 2D) 



at frequency regime lo > (<)(1/to, fc^T/S), where tq is 
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the elastic scattering time 18| . The dielectric function is 
given at small q, uj by 

£(q, 0.) = 1 - ^Xd(q, c^) ^ 1 + 47r/3e2 + 0{q^), (7) 

q 

also in agreement with the result obtained in U(l) gauge 
theory [l8|- 

The U(l) spin liquid is actually a member of the QSLs 
described by our phenomenology. To show this we start 
with a Landau Fermi liquid with interaction parameters 
Fq and F^{q,u!) only. The long- wavelength and low dy- 
namics of the Fermi liquid is described by an effective 
Lagrangian 



^eff ■ 



E 



(8) 



where c^g.(ck(T) are spin-cr fermion creation (annihilation) 
operators with momentum k, and 



^ ' ' 2N{0) ^ 



(9) 

describes the current-current and density-density inter- 
actions between quasi-particles [l7l| . where q = (q, w) and 
vp = hkp/m* is the Fermi velocity. 

The current- and density- interactions can be decou- 
pled by introducing fictitious gauge potentials a and if 
(Hubbard-Stratonovich transformation) with 
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J • a + n^ - - ——a 

2 \m* Fl 
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N{0) 2 
F^ ^ 



(10) 

where n is fermion density. We have used the equality 
d{n/m*) = N{0)vp in writing down Eq. pH)) . 

The Lagrangian ([5]) and ([TU]) can be rewritten in the 
standard form of U{1) gauge theory by noting that the 
fermion current is given in this representation by 



2m' 



where ipai"^) — /e ^ 'cko- in the Fourier transform of 
Ckcr- The Lagrangian can be written as 



where 
and 



d'^r 



^t(j^_^)^^-i/(^t,V.) 



H(V't,V.) = :— ^|(V-^a)^.p 
2m* 



+ L{ip, a), 
(11a) 

(lib) 



L((p,a) = i J d'l 



■^(1 + A)a2 
m*^ F{' 



Fri ^ 



(lie) 



Notice how the 5^^^ term in L{(p,BL) arises from the 
introduction of diamagnetic term in H{ipl., ipa)- 

Using Eq. ([5|), we find that in the small q limit, the 
transverse part of L{ip, a) is given in the spin liquid state 
by 



2m* 



d"" 



7t(V X a)2 



(12) 



The longitudinal part of the gauge potential ip is screened 
as long as Fq ^ 0. Lagrangian PT|) together with ([T^ is 
the standard Lagrangian used to describe C/(l) QSLs. It 
is interesting to note that a nonzero 1 + ^^i^ leads to a 
mass term for the gauge field a, in agreement with slave- 
boson/rotor approaches where a metallic state appears 
with condensation of bosons/rotors [I6, 19]. 

Our picture of QSLs has several immediate experimen- 
tal consequences. For example, the specific heat ratio 
7* — Cv/T = m*/m.Cy^ is predicted to be finite at 
the QSLs where Cy^ is the specific heat for the corre- 
sponding non-interacting electron liquid. The magnetic 
susceptibility given by 



1 



XP 



_ (0) 



m 1-f Fq" 

is also non-singular in the QSLs with Wilson ratio R = 

Zlfrr (1 + ^Sr' generally of order 0(1). The 
transport properties in the QSLs can be computed us- 
ing the Landau transport equation. We expect that the 
transport life-times will be dominated by scattering in 
the current-current channel which is the most singular 
scattering channel in the limit 1 + -j- 0. The effect 
of scattering in this channel should be equivalent to re- 
sults obtained from U{1) gauge theory. The thermal con- 
ductivity K computed from Landau transport equation is 
given by 



K 

T 



1 



TVp To 



where I/tq is the elastic (impurity) scattering rate, ep — 
p'^/2m* is the spinon Fermi energy and vp is the spinon 
Fermi velocity. The same result is obtained in U(l) gauge 
theory at 2D.[23l (See supplementary material for de- 
tails.) 

The close relation between Fermi liquid and spin- 
liquid states suggests that the (zero temperature) metal- 
insulator transition between the two states is character- 
ized by the change of Landau parameter 1 + Ff{0, 0)/d — > 
0''" across the transition. The nature of metal-insulator 
transition within the Fermi liquid framework was first ad- 
dressed by Brinkman and Rice[2lj where they proposed 
that a metal-insulator (Mott) transition is indicated by 
diverging effective mass ^ >■ 00 and inverse compress- 
ibility K — at the Mott transition point, with cor- 
respondingly a vanishing quasi-particle renormalization 
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weight Z ~ ^ 0. The diverging effective mass and 
vanishing quasi-particle weight suggest that the Fermi 
hquid state is destroyed at the Mott transition, and the 
Mott insulator state is distinct from the Fermi hquid 
state at the metal side. Here we propose an alterna- 
tive picture where the Fermi surface is not destroyed but 
the quasiparticles are converted into spinous at the Mott 
transition. In particular, the effective mass m*/m may 
not diverge at the metal-insulator transition although 
Z — > in this picture. A schematic phase diagram for 
the Mott (metal-QSLs) transition is presented in Fig. ([T]) 
where we imagine a Hubbard type Hamiltonian with hop- 
ping t and on-site Coulomb repulsion [/. The system is 
driven to a Mott insulator state at zero temperature at 
U = [/c, where 1 -f- F{{U > Uc)/d = 0. We caution 
that in general a finite (T 7^ 0) region exists around the 
Mott transition point where the physics is dominated by 
critical fluctuations and our phenomenological theory is 
not applicable. We note that an alternative phenomenol- 
ogy for Mott transition from a semi-microscopic start- 
ing point [22] has qualitatively similar conclusion as our 
present work. The relation between the two works is not 
clear at present. 
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FIG. 1: (a) Schematic zero temperature phase diagram 
for Mott transition. U is the Hubbard interaction strength 
and t is tlie hopping integral. The electron quasiparticle 
weight and quasi-particle charge current ~ 1 -I- Ff /d van- 
ishes at the critical point while the effective mass remains 
finite, (b) Schematic phase diagram showing finite tempera- 
ture crossovers and possible instability toward gapped phases 
at lower temperature. There exists a (finite temperature) crit- 
ical region around Uc where our phenomenological theory is 
not applicable. 



We end with a comment about the Pomeranchuk in- 
stability. Experienced researchers in Fermi liquid theory 
will recognize that the point 1 + F^/d = is in fact a 
critical point in Fermi liquid theory. The Fermi surface is 
unstable with respect to deformation when 1 + Ff/d < 0. 
The stability of the 1 + Ff/rf = point is required in 
QSLs where quasi-particles (spinous) become chargeless. 
The resulting QSLs we obtain here are marginally stable 
because of large critical fluctuations. Similar result was 
obtained in U{1) gauge theory, where it was found that 
a fermionic spin liquid state with large Fermi surface is 
marginally stable toward C/(l) gauge fluctuations. The 
Pomeranchuk criticality is an alternative way to express 
this result. 

The presence of Pomeranchuk criticality suggests that 
QSLs with large Fermi surfaces are in general rather sus- 
ceptible to formation of other more stable QSLs at lower 
temperature, like the Z2 QSLs or valence bond solid 
(VBS) states that gap out part or the whole Fermi sur- 
face. The resulting phase diagram at the vicinity of Mott 
transition thus has the generic feature shown in Fig.([TlD), 
where the system is driven into a gapped QSL at low tem- 
perature T < Fc{U) at the insulating side. The nature 
of the low temperature QSLs depends on the microscopic 
details of the system and cannot be determined from our 
phenomenology. Our theory is applicable at T > Fc{U), 
when the spin liquid is still in the large-Fermi surface 
phase. 
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RENORMALIZED CURRENTS 

In this supplementary material we derive the renor- 

malizcd currents in Fermi liquid theory ((Eq.(2)) in the 
main text). The local equilibrium quasi-particle occupa- 
tion numbers and their fluctuations have to be considered 
carefully. The excitation energy of an additional quasi- 
particle with momentum p is given by 



local equilibirum), 



= ep + X]-^PP'^"p'' 



where ep ~ . The corresponding local equilibrium oc- 
cupation number is hp = np (cp — A*), and the departure 
from local equilibrium reads 



Sh„ 



p' 



where Sup = rip — rip. The charge current J carried by 
quasiparticles is related to the particle density by the 
conservation law, 



0. 



The density fluctuation Sp (r, t) should be expressed in 
terms of the sum of Sfip (r, t) (i.e. fluctuation away from 



and 



Therefore, 



where 



—6p + Vr • ^ ShpVp = 0. 



Jp = Vp-2^/pp'^Vp.. 



p' 



Using the relation 



Err/pp'Vp' = J2^^f'J:piP^(^^'^)^^ 



p' 



p' 



1 "pi 



where Np = N (0), we flnd that the renormalized charge 
current is 

J = ^(l + iFf)jW, 

when J*^ is the electric current in the absence of inter- 
action. Notice that Pi (cos 6) is replaced by cos {W) at 
2D. 

Similarly the thermal (energy) current Jg is given by 

= XI ^"p *^^p ~ '")^p 

p 

= X] (^P ~ ^) Vp(5np - X |r-/pp"5"p') 



p' 



X] 5np[{ep - m) Vp - X] 1^7 -^pp' (^p' ~ ^p']- 



p' 



Notice that 



XIItT^pp' (v -M)vp' 
p' ^ 

X Ti^^F ' E ^iP^ (^°« ^) (^p' - ^p' 

p' I 



1 f dn" 



to leading order. Therefore the renormalized energy cur- 
rent is given by 

We observe that the thermal current is not renormalized 
by the factor (l + iff). 
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LANDAU TRANSPORT EQUATION 



Using the expansion 



In the Landau transport equation which is essen- 
tially a Boltzmann equation, the transition probability 
W{1, 2; 3, 4) for a two quasi-particles scattering process in 
an isotropic Fermi liquid, 1 + 2—^3-1-4 with i = (p^, ct^), 
is given by 27r times the squared moduli of the quasi- 
particle scattering amplitude. We are interested in the 
situation that the momentum transfer q = pi — p3 is 
small and p = ^(pi + Ps) and p' — i(p2 + P4) are close 
to the Fermi momentum p^^. In this case, the scattering 
amplitude depends mainly on the relative orientation of 
the vector p, p' and q, and on the the energy transfer 

Quasi-particle scattering amplitude 

We shall neglect the spin indices in the following for 
brevity. The spin-degeneracy factor 2 will be inserted 
when the need arises. The quasi-particle scattering am- 
plitude Appi (q, Lu) is then given by the following equa- 
tion, 

App' (q,cj) - ^/pp"Xop" (q,w)-4p"p' (q,w) = /pp', 
p" 

where xop (l: ^) is the susceptibility 



xop (q, t^) = 



'■p-q/2 



,0 

''p+q/2 



q ■ Vp ^'^p 

Cp-q/2 - Q ' Vp - W 9^p ' 



with = ni?(^k)- We shall assume that the scattering 
is dominating by the I = 1 channel and approximate 



JPP' „2 ^1- 

It is then easy to show that 
P P' 



^PP' 



ft 



p% i-xi(q,w)/f 



where 



Xi 



_ ^ 

"-p-q/2 ~ "-p+q/2 
'^ + Cp-q/2 



For small (q, a;) we have 
Xi(q,w) 



NiO) 



l + tg{d) 



for LU ^ vpq, where q = |q|, g{2) = 1 and g{3) ~ ^. 
In the limit N{0)fl/d = Ff/d -1, the scattering 
amplitude becomes 



App, (q, w) 



P P 



Pf l-Xi(q,w)/f 



we obtain 

^pp' (q.w) 



N{0) 

d pp' 1 

N{0) pI -z.g^+7,g2' 



(13) 



where the last factor is exact the gauge field propagator 
in U{1) gauge theory. 



Thermal conductivity 

Following Pethick[l], we shall use a variational 
approach [2] to derive the thermal conductivity k from the 
transport equation. The thermal resistivity for a Fermi 
liquid is given by 

- = I E W^(l,2;3,4Kn°(l-n°)(l-n^) 



1,2,3,4 



X {(f)i -1-02-03-0' 



X(5(ei + £2 - £3 - e4)(5CTi+(T2,cr3+CT4(5pi+p2,P3+pi>l-4) 

where n? — np{^i) is the Fermi distribution function with 
i = 1,2,3,4, = ti — /i, 4>i is defined by rii = — 

4>i-Q^, Vi is the quasiparticle velocity, u is an arbitrary 
unit vector along which the temperature gradient VT is 
applied. 

We shall take the usual trial function = ^^Vi ■ u. To 
the order we are working on the quasi-particle velocity 
may be replaced by and 



4m*iV(0) Z"^^,^ , ..M2dn^{0 



4TO*iV(0) TT 



27r n 
3 m* 



{kBT)\ 



where n is the fermion density and the relation 
d{n/m*) = N{0)vl is used. 

Introducing ^p = ^{£,p+q/2 + Cp-q/2) and using the 
conditions q = pi — P3 = P4 — P2 and u = ep_|_q/2 — 
ep~q/2 ep,+q/2 - ep'-q/2, wc have 

(ClVl + 6V2 - C3V3 - C4V4) 

= (Cp + oj/2) (p + q/2) + - Lo/2) (p' - q/2) 
- (Cp - C./2) (p - q/2) - (fp, + C./2) (p' + q/2) 

= (Cp-Cp')q+'^(p-p') 
= (^p -Cp')q+w(p-p')> 
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and 



{{4>1 + (j)2 - 04) ) 

([(^Vi + ^2V2 - - ?4V4) • U]' 

1 2 
^ (^iVl + C2V2 - - ^4V4) 

^ [(ep-^p')q+(p-p'H', 



where (• • • ) means averaging over different u. 

Putting App> {q,uj) into P^ . and using the identity 
nO(l - n°) = {n° - - e^^'^-^^^, we obtain 



Assuming that the scattering is dominated by channel 
and using Ea. (jl3p . we obtain 



1 



dio- 



1 

7 "'^(e'3"-l)(l-e-'3-')^ 



-(4-d)/3 



1 1 

— cx 



y4 5Z 

q.p.p 



I^PP' (q,^)l' 



- l)(l e-'3") 

x[nF(Cp-q/2) - "F(^p+q/2)]'5(w - P ' q/w*) 
x["F(Cp'+q/2) - "F(Cp'-q/2)]'5(w - p' ■ q/m*) 



x[(Cp-^P')V+(p-p'M, 



where we have used the (5-functions to simphfy the 
last line. Replacing nF(^p-q/2) - '^F(Cp+q/2) and 
'^F(?p'+q/2) - «F(Cp'-q/2) by and respec- 

tively which is valid at small q and w, we obtain 



which is the thermal resistivity coming from inelastic 
scattering between fermions. At low temperature the in- 
elastic scattering is cut off by elastic impurity scattering 
rate I/tq which gives rise to 



1 



where 7* — Cy/T is the specific heat ratio. The total 
thermal conductivity is therefore given by 



q,p,p' 



UJ^\App, (q,w) p 



(e/^-^ - 1)(1 -e-/3'^) 
S{u} — p • c[/m*)S{uj — p' • q/™*) 



x[(ep-ep')V+(P-p'M- 

Let 0(6') be the angle between p(p') and q, and integrat- 
ing over and ^p' , we obtain 



1 1 
— oc 

K 



rp4 $Z 

q,p,p 



{el^'^ - 1)(1 -e-^") 
x6{bj — qvp cos9)5{ijj — qvp cos9') 



K 

— oc max 
T 



fr2 



EF 



1 



T^F To 
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